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ABSTRACT 

We use time domain analysis techniques to investigate the rapid variability of Cygnus X-l. We show that the 
cross-correlation functions between hard and soft energy bands reach values very close to unity and peak at a lag of 
less than 2 millisecond for energies separated by a factor of 10. This confirms that the process that produces X-ray 
photons at different energies is extremely coherent on short time scales and strongly constrains emission models 
proposed to explain Fourier-frequency-dependent time lags. We present autocorrelation functions at different en- 
ergies, and note their widths decrease with increasing energy. We show that the extended Compton corona model 
produces auto-correlation functions whose widths increase with increasing energy, that the model of cylindrical 
waves moving inward through a transition disk has too large a peak lag in the cross-correlation function. Models 
of magnetic flaring and of drifting blobs in a hot corona can qualitatively fit the observations. 

Subject headings: accretion, accretion disks - black hole physics - methods: data analysis - stars: individual 
(Cygnus X-l) - X-rays: stars 
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1. INTRODUCTION 

The X-ray/7-ray spectrum of an accreting black hole such as 
Cygnus X-l in its hard state can be represented as the sum of 
a few components: a soft component associated with the emis- 
sion from a cold accretion disk, a hard tail extending up to a 
few hundred keV associated with a hot "corona", and a Comp- 
ton reflection bump produced when hard X-rays are reflected 
from cold material in the accretion disk (see, e.g., Zdziarski 
et al. 1997; Gierlinski et al. 1997; Poutanen 1998). Spectral 
data suggest thermal Comptonization by a medium with a tem- 
perature of about 100 keV as the origin of the hard tail. The 
observed X-ray spectral slopes and the amplitude of Compton 
reflection can be used to determine the geometry of the sys- 
tem. However, a variety of models fit the spectra (see Poutanen 
1998; Beloborodov 1999a; Zdziarski 2000), so the parameters 
of the accretion flow are not be well constrained by spectral 
data alone. Most of the spectral models, however, are applied to 
time-averaged spectra despite of the fact that the sources show 
rapid spectral variability (see, e.g., Nolan et al. 1981; Negoro, 
Miyamoto, & Kitamoto 1994; Feng, Li, & Chen 1999) imply- 
ing rapid changes in the physical conditions. Examining tem- 
poral characteristics should help break the degeneracy among 
spectral models. 

Time domain techniques were used in the early days of X- 
ray astronomy before the statistical samples of data were suffi- 
cient to make use of Fourier domain analyses (see, e.g., Weis- 
skopf, Kahn, & Sutherland 1975; Sutherland, Weisskopf, & 
Kahn 1978; Pnedhorsky et al. 1979; Nolan et al. 1981). The 
asymmetry of the cross-correlation function (CCF) of Cygnus 
X-l was discovered in ~ 150 sec observations by Priedhorsky 
et al. (1979) and Nolan et al. (1981). They showed that the CCF 
peaks at a lag < 10 — 40 ms. Using data from EXOSAT, Page 
(1985) confirmed these results and claimed a ~ 6 ms shift of the 
peak of the CCF between the 5-14 keV and the 2-5 keV bands. 
These are the last papers, to our knowledge, that present the 
CCFs despite immense advances in temporal resolution, photon 



statistics, and duration of observations. Aside from attempts to 
model individual shots (Lochner, Swank, & Szymkowiak 1991; 
Negoro et al. 1994; Focke 1998; Feng et al. 1999), recent anal- 
yses have concentrated on Fourier domain techniques. 

The CCF asymmetry is related to the Fourier-frequency- 
dependent hard time lags between different spectral bands dis- 
covered by Ginga (Miyamoto et al. 1988). These data gave new 
strong constraints on spectral models. However, in the Fourier 
domain it is difficult to measure time lags at frequencies above 
~ 30 Hz (see, e.g., Nowak et al. 1999a) corresponding to the 
light travel time in a region of the main energy dissipation. On 
the other hand, in the time domain, the CCFs can be measured 
accurately down to the lags as short as 2 ms. Furthermore, time 
domain functions (containing, in principle, the same informa- 
tion as their Fourier domain companions) highlight different 
information which can be used to constrain emission models 
further. 

We present the results of time domain analyses of Cygnus X- 
1 as observed by the Rossi X-Ray Timing Explorer (RXTE) and 
compare our observational results with model predictions. 

2. OBSERVATIONS AND RESULTS 

2.1. Observation Log 

Cygnus X-l was observed in its hard (low) state 12 times by 
RXTE during December of 1997, for a total of about 30 ksec. 
Here we present results only from the Proportional Counter 
Array, using the standard screening criteria of earth elevation 
greater than 10°, offset from source less than 0.01°, all 5 pro- 
portional counter units on, and the standard time since the 
last South Atlantic Anomaly passage. For the lowest ener- 
gies (below 8 keV), we have data recorded in single bit (SB) 
modes, where RXTE counts photons with no spectral infor- 
mation other than whether they fall within the given range of 
channels. For the higher photon energies, we use the full Event 
Mode data with all the available spectral information. The use 
of SB modes and the subsequent loss of some spectral informa- 
tion were required because of telemetry limitations of RXTE. 
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FIG. 1 . — The auto- and cross-correlation functions of Cygnus X-l observed in December 1997. Solid curves show the ACFs for the 2-5 keV energy band (thick 
curves) and the 24-40 keV band (thin curves). Dotted curves show the CCF between the 8-13 keV band and the 2-5 keV band, and dashed curves represent the CCF 
for the 24-40 keV band vs the 2-5 keV band. The CCFs at other energies have similar shapes. The peaks all align at around zero lag. The higher energy curves are 
narrower. The CCFs are defined in such a way that the peak is expected to appear at a positive lag when hard photons are lagging the soft ones. The errors at small 
lags are roughly the same size as the line widths and are hence left unplotted for clarity. The CCF peak lag is at less than 2 ms for the 2-5 keV vs 8-13 keV CCF. 



We make no background subtractions for our data because the 
background counts are a very small fraction of the total counts 
and because current background models for RXTE do not give 
estimates on time scales shorter than 16 seconds. RXTE dead- 
time for Cygnus X-l is about 1 % and should affect only zero 
lag bins in the correlation functions. 

2.2. Light Curves 

Using the standard FTOOLS 4.1 software, we extract light 
curves with 2 -9 second resolution (~ 2 ms) for each of the SB 
modes (for energy channels 2-5 keV, 5-6 keV, and 6-8 keV), 
plus the event mode data binned into three additional energy 
bands (channels 24-40, 41-71, and 72-132, or 8-13 keV, 13- 
24 keV, and 24-40 keV respectively). We do not analyze any 
higher energy data because the background count rates become 
large above 40 ke V. We attempted to analyze the data with a ~ 1 
ms bin size. For the low energy (and hence higher count rate) 
bands the 1 ms computations allowed us to produce slightly 
stronger constraints (discussed below). We plot only the 2 ms 
binning results in order to make the plots clearer, since the re- 
sults are essentially the same in either case. 

2.3. Cross-correlation Functions 

We compute the CCF of Cygnus X-l, comparing each energy 
band with the lowest energy band. The CCFs are asymmetric 
with peaks at lags less than 2 ms in all cases (see the inset of 
Fig. 1), and less than 1 ms between the 3 keV band and bands 
below 13 keV. The rising part of the CCF (soft lags) becomes 
narrower with energy substantially faster than the decaying part 
(hard lags). These results agree qualitatively with past results 
for the CCF (Priedhorsky et al. 1979; Nolan et al. 1981). The 
CCFs reach values very close to unity, showing that the signal 
at all the energies is extremely well synchronized. The asym- 
metry of the CCFs in the time domain has a direct relation to 



the Fourier- frequency-dependent time lags. 

2.4. Autocorrelation Functions 

The autocorrelation functions (ACFs) of Cygnus X-l are 
shown in Figure 1 and Figure 2a. The width of the ACF de- 
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creases with photon energy approximately as oc E 
lags smaller than ~ 0.3 sec. (At larger lags the ACFs at dif- 
ferent energies are not self-similar.) This strongly constrains 
the origin of the spectral variability, since it requires that the 
pulses producing the variability last longer at low energies than 
at higher energies. Similar energy dependence is observed also 
in the ACF of the peak aligned average shot profiles (Feng et 
al. 1999). Our results extend their work by demonstrating the 
trend of width of the ACF versus energy applies across more 
energy bands. More importantly, we prove that this trend is 
always present in the data, and is not subject to the selection 
effects of a shot-fitting algorithm. 

3. CONSTRAINTS ON THE EMISSION MODELS 

3.1. Implications on Shot Shapes for Shot Noise Models 

One can attempt to explain the observed ACFs and CCFs in 
terms of a simple shot noise model (Terrell 1972). Since the 
CCFs peak at a lag < 2 ms, the shots at different energies also 
should reach maxima within 2 ms from each other. We can 
develop some intuition about the typical shot profiles by look- 
ing at the analytic form of the CCF for shots which have expo- 
nential rise and decay and peak at the same time for different 
energies. The CCF is 

r s 2 d (ri lr + r hd )f7/^ T h 2 r (r sr +r sd )e^ 



oc 



(Thd + T s d)(T s d - Thr) (Thr + T sr )(T sd - T hr ) 

for negative r, and 

(r sr + T sd )e- T /^ d T s 2 r (r hr + r hd ) e - T / T " 
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Fig. 2. — The ACFs in the different energy bands, (a) The data for Cygnus X-l. Solid curves shows the ACF for the 2-5 keV energy band, dotted curve is for 8-13 
keV band, and dashed for 24-40 keV band, (b) The prediction of the extended corona model. 



for positive r, where T sr ,T s d and Th r ,Thd are the rise and 
decay time constants of the shots in soft and hard en- 
ergy bands, respectively. The proportionality coefficient is 
2/ sj (r sr + T s d)(Thr + Thd)- The energy dependence of the 
ACF requires the shots at higher energies to be shorter, i.e. 
Thd + "Tir < 7sr + T"sd- Since the observed CCF (see Fig. 1) is 
roughly equal to the ACF of the softer energy band for positive 
r and the ACF of the harder energy band for negative t, and 
the CCF rises faster than it decays, one gets max(r s d, Th r ) ~ 
max(Th r ,Thd) < max(r sr ,Thd) ~ max(r sr , r s d). This means 
that the soft rise time scale is the longest, r sr > r s d, T[ lr , Thd, 
and that the decay times at different energies are very close to 
one another (i.e., T s d ~ Thd) or are so small that the only rele- 
vant parameters are the rise times (Th r 3> T s d, Thd). The energy 
dependence of the rise time then produces the time lags and the 
asymmetry of the CCF (see also Miyamoto & Kitamoto 1989 
who arrived at similar conclusions). 

In modified shot noise models (see, e.g., Lochner et al. 
1991), there is a broad distribution of shot time scales. The 
shots at different energies should be perfectly synchronized in 
order to achieve high values of the CCFs. However, constraints 
on the shape of the shots are not so strong as for the simplest 
shot noise model. The shots at different energies can have simi- 
lar properties as described above, or they can be shifted relative 
to each other (e.g., Poutanen & Fabian 1999a,b). If the shift de- 
pends on the shot time scale as oc r a , one obtains Fourier time 
lags St(f) cx f~ a . In such a case, the CCF peaks at a lag equal 
to the delay corresponding to the shortest time scale. The data 
then constrain the minimum shot time scale to be < 1 ms. 

3.2. Extended Corona Models 

Comptonization in a uniform electron cloud produces time 
lags which are frequency independent in the range we can probe 
with current X-Ray instruments (Miyamoto et al. 1988). This 
inspired Kazanas, Hua, & Titarchuk (1997) to propose an ex- 
tended corona model with a r _1 radial density distribution and 
a size of a few light seconds. The long lags observed at low 
Fourier frequencies are produced here by photons travelling and 



scattering over large radii, while the shorter time lags are pro- 
duced in the central small core of the cloud. The model fits 
much of the data, but has the physical problem of having too 
much energy input at large radii. 

We simulated the light curves for this model and computed 
the auto-correlation functions at different energies. As an ex- 
ample, we consider a set of parameters from Hua, Kazanas, 
& Cui (1999) (kT r = 0.2 keV, kT e = 100 keV, p = 1, 
m = m = 10 16 cm- 3 , n = 10~ 3 lt-s, r 2 = 10 3 ri, for 
details see discussion around their Fig. 1). The shots are as- 
sumed to be produced by modulations in the soft flux. Because 
this model produces delays by having larger light travel times 
for higher energy photons, it always produces a wider shot and 
hence a wider ACF at higher energies (see Fig. 2b). We note 
that the ACF at 30 keV in this model is a factor of ten broader 
than the ACF at 3 keV, while in the observations it is 50% nar- 
rower. The fact that the ACFs become broader with energy is 
the intrinsic property of the model and there is no way to re- 
solve this problem by changing the parameters of the system. 
We conclude that models where the time lags are produced by 
light travel delays can be ruled out. 

3.3. Magnetic Flare Models 

Magnetic flares on the surface of the cold accretion disk were 
shown recently to produce XAy-ray spectra in agreement with 
observations of Cyg X-l (Beloborodov 1999b). The observed 
time lags can correspond to the time scale of the evolution of 
magnetic structures (Poutanen & Fabian 1999a,b). The mag- 
netic field lines twist due to differential rotation in the accretion 
disk and elevate to the corona releasing magnetic energy and 
heating the corona. The flare time scale is then of the order of 
the Keplerian time scale at the relevant distance from the cen- 
tral black hole. Changes in the energy dissipation rate and in 
the geometry of the flare (e.g., distance from the disk) produce 
soft-to-hard spectral evolution which is the cause of the hard 
time lags. A broad distribution of the flare time scale assures 
that the time lags are inversely proportional to the Fourier fre- 
quency St(f) oc Tf ~ l/(27r/), where t/ is the shot time scale 
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giving contribution to the power spectrum at frequency /. 

The model of Poutanen & Fabian (1999b) reproduces well 
the time lags observed in Cyg X-l, while producing a somewhat 
wider ACF at higher energies due to the assumption that the en- 
ergy dissipation rate rises faster than it decays. In the opposite 
situation, when the dissipation rate rises more slowly than it de- 
cays, the ACF and CCF energy dependences can be reproduced 
easily (see, e.g., Miyamoto & Kitamoto 1989; Poutanen 2000). 

3.4. Drifting Blob Models 

Bottcher & Liang (1999) proposed a model in which spectral 
variability is produced by a cool blob drifting inward through 
an inhomogeneous hot inner disk. This model qualitatively 
matches the data. However, the parameters presented in that 
work, produce significant quantitative deviations from the ob- 
servations in both the Fourier and time domains. While chang- 
ing the parameters of the system could allow the model to fit 
the data, a larger problem for this model is that it drives the 
variability through modulations of the soft photon flux. If the 
energy dissipation rate in the corona does not change with time, 
then an increase of the soft photon flux would lead to a softer 
spectrum and spectral pivoting around ~ 10 keV. The ampli- 
tude of the variability would be then a strong function of photon 
energy (larger at larger energies) and the variability above and 
below the pivoting point would anti-correlate, contrary to what 
is observed (see, e.g., Nowak et al. 1999a and § 2.3). In or- 
der to reproduce the observed correlated variability at different 
energies without violating the energy balance, one has to as- 
sume that the inward drift of cool blobs is perfectly correlated 
with the increase of the energy dissipation in the hot corona. It 
remains to be seen whether such a requirement is physically re- 
alistic and the model can indeed fit the data. One cannot avoid 
this constraint by proposing an ADAF-type solution where the 
coronal cooling rate is dominated by cyclo-/synchrotron radi- 
ation rather than by Compton cooling. Such a model would 
require a cyclo-/synchrotron luminosity (observed in IR/optical 
light) far in excess of the X-Ray luminosity of ~ 10 37 ergs/sec. 
The observed optical luminosity is only about 5 x 10 36 ergs/sec 
and is dominated by the companion star. 

3.5. Cylindrical Wave Models 

Cylindrical waves propagating through the accretion disk 
from the region where the soft X-rays are emitted to the re- 



gion where the hard X-rays are emitted (e.g., Miyamoto et al. 
1988; Kato 1989; Nowak et al. 1999b) have also been proposed 
as a mechanism for producing hard time lags. The dispersion 
of wave velocities results in the frequency dependence of the 
time lags. A recent and relatively well developed version of 
such a model (a "transition disk" model) designed to explained 
the spectrum of Cyg X-l and the Fourier time lags was recently 
considered by Misra (2000). However, in this model the peak 
of the emission at 30 keV is delayed by ~ 0.015 s from the 
peak at 3 keV. In such a case, the CCF between the 30 keV and 
3 keV photons would peak at that lag strongly contradicting the 
data. One can assume that the low frequency signal propagates 
slowly producing larger time lags, while high frequency signal 
propagates faster producing smaller time lags. (This is basically 
a modified shot noise model.) The transition disk model then 
requires a propagation speed of > 2c in order to fit < 2 ms lag 
in the CCF. 

4. CONCLUSIONS 

We present the results of cross-correlation and autocorrela- 
tion analysis of the light curve of Cygnus X-l in several energy 
bands. The width of the ACF was found to become narrower 
with photon energy as oc E~°- 2 . The corresponding CCFs are 
asymmetric, but all peak at lags less than 2 ms. 

We compare these results to model calculations for an ex- 
tended Compton corona model, a magnetic flare model, a drift- 
ing blob model, and a cylindrical wave model. The extended 
corona model inherently produces longer shots at higher ener- 
gies (and, therefore, broader autocorrelation function), and it 
cannot fit the data even qualitatively. We find that a magnetic 
flare model can fit all the data if one requires that the energy 
dissipation rate rises slower than it decays. We also find that 
while the drifting blob model fits the data qualitatively, it re- 
quires the inward drift of cool blobs to be perfectly correlated 
with the increase of the energy dissipation in the hot corona. A 
transition disk model predicts a large shift of the peak of the 
cross-correlation function from zero lag contradicting the data. 

We thank Steve Kahn for useful discussion and Andrei Be- 
loborodov for valuable comments. This research has been sup- 
ported by NASA grant # NAGS-6691 and NAGS-7409 (TJM, 
PSC), and by the Swedish Natural Science Research Council 
and the Anna-Greta and Holger Crafoord Fund (JP). 
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